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25.1 Parametrizing the Universe 

Rapid advances in observational cosmology have led to the establishment of a precision cosmo- 
logical model, with many of the key cosmological parameters determined to one or two significant 
figure accuracy. Particularly prominent are measurements of cosmic microwave background (CMB) 
anisotropies, with the highest precision observations being those of the Planck Satellite [1,2] which 
supersede the landmark WMAP results [3,4]. However the most accurate model of the Universe 
requires consideration of a range of observations, with complementary probes providing consistency 
checks, lifting parameter degeneracies, and enabling the strongest constraints to be placed. 

The term ‘cosmological parameters’ now has a wide scope, and may include the parameterization 
of some functions as well as simple numbers describing properties of the Universe. The original usage 
referred to the parameters describing the global dynamics of the Universe, such as its expansion 
rate and curvature. Now we wish to know how the matter budget of the Universe is built up 
from its constituents: baryons, photons, neutrinos, dark matter, and dark energy. We also need 
to describe the nature of perturbations in the Universe, through global statistical descriptors such 
as the matter and radiation power spectra. There may be additional parameters describing the 
physical state of the Universe, such as the ionization fraction as a function of time during the 
era since recombination. Typical comparisons of cosmological models with observational data now 
feature between five and ten parameters. 


25.1.1 The global description of the Universe 

Ordinarily, the Universe is taken to be a perturbed Robertson-Walker space-time, with dy- 
namics governed by Einstein's equations. This is described in detail in the Big-Bang Cosmology 
chapter in this volume. Using the density parameters Q; for the various matter species and (2, for 
the cosmological constant, the Friedmann equation can be written 
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where the sum is over all the different species of material in the Universe. This equation applies at 
any epoch, but later in this article we will use the symbols Q; and Q, to refer specifically to the 
present-epoch values. 

'The complete present-epoch state of the homogeneous Universe can be described by giving the 
current-epoch values of all the density parameters and the Hubble constant h (the present-day 
Hubble parameter being written Hy = 1005 km s^! Mpc ^ !). A typical collection would be baryons 
Qp, photons Q}, neutrinos Qy, and cold dark matter Q. (given charge neutrality, the electron 
density is guaranteed to be too small to be worth considering separately and is effectively included 
with the baryons). The spatial curvature can then be determined from the other parameters using 
Eq. (25.1). The total present matter density Qm = Ne + Np may be used in place of the cold dark 
matter density Qe. 

'These parameters also allow us to track the history of the Universe, at least back until an epoch 
where interactions allow interchanges between the densities of the different species; this is believed 
to have last happened at neutrino decoupling, shortly before Big-Bang Nucleosynthesis (BBN). To 
probe further back into the Universe's history requires assumptions about particle interactions, and 
perhaps about the nature of physical laws themselves. 
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The standard neutrino sector has three flavors. For neutrinos of mass in the range 5 x 1074 eV 
to 1 MeV, the density parameter in neutrinos is predicted to be 


2 My 

Oh? = A LT (25.2) 
where the sum is over all families with mass in that range (higher masses need a more sophisticated 
calculation). We use units with c = 1 throughout. Results on atmospheric and Solar neutrino 
oscillations [5] imply non-zero mass-squared differences between the three neutrino flavors. These 
oscillation experiments cannot tell us the absolute neutrino masses, but within the normal assump- 
tion of a mass hierarchy suggest a lower limit of approximately 0.06eV for the sum of the neutrino 
masses (see the Neutrino chapter). 

Even a mass this small has a potentially observable effect on the formation of structure, as 
neutrino free-streaming damps the growth of perturbations. Analyses commonly now either assume 
a neutrino mass sum fixed at this lower limit, or allow the neutrino mass sum to be a variable 
parameter. To date there is no decisive evidence of any effects from either neutrino masses or 
an otherwise non-standard neutrino sector, and observations impose quite stringent limits; see the 
Neutrinos in Cosmology chapter. However, we note that the inclusion of the neutrino mass sum as 
a free parameter can affect the derived values of other cosmological parameters. 


25.1.2 Inflation and perturbations 

A complete model of the Universe must include a description of deviations from homogeneity, 
at least statistically. Indeed, the most powerful probes of the parameters described above come 
from the evolution of perturbations, so their study is naturally intertwined with the determination 
of cosmological parameters. 

There are many different notations used to describe the perturbations, both in terms of the 
quantity used to and the definition of the statistical measure. We use the dimensionless power 
spectrum A? as defined in the Big Bang Cosmology section (also denoted P in some of the lit- 
erature). If the perturbations obey Gaussian statistics, the power spectrum provides a complete 
description of their properties. 

From a theoretical perspective, a useful quantity to describe the perturbations is the curvature 
perturbation R, which measures the spatial curvature of a comoving slicing of the space-time. 
A simple case is the Harrison—Zeldovich spectrum, which corresponds to a constant AS. More 
generally, one can approximate the spectrum by a power law, writing 


k ng—1l 
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where ns is known as the spectral index, always defined so that ns = 1 for the Harrison—Zeldovich 

spectrum, and k, is an arbitrarily chosen scale. The initial spectrum, defined at some early epoch 

of the Universe's history, is usually taken to have a simple form such as this power law, and we will 

see that observations require ns close to one. Subsequent evolution will modify the spectrum from 

its initial form. 

'The simplest mechanism for generating the observed perturbations is the inflationary cosmology, 
which posits a period of accelerated expansion in the Universe's early stages [6, 7]. It is a useful 
working hypothesis that this is the sole mechanism for generating perturbations, and it may further 
be assumed to be the simplest class of inflationary model, where the dynamics are equivalent to 
that of a single scalar field ¢ with canonical kinetic energy slowly rolling on a potential V(¢). One 
may seek to verify that this simple picture can match observations and to determine the properties 
of V(@) from the observational data. Alternatively, more complicated models, perhaps motivated 
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by contemporary fundamental physics ideas, may be tested on a model-by-model basis (see more 
in the Inflation chapter in this volume). 

Inflation generates perturbations through the amplification of quantum fluctuations, which 
are stretched to astrophysical scales by the rapid expansion. The simplest models generate two 
types, density perturbations that come from fluctuations in the scalar field and its corresponding 
scalar metric perturbation, and gravitational waves that are tensor metric fluctuations. The former 
experience gravitational instability and lead to structure formation, while the latter can influence 
the CMB anisotropies. Defining slow-roll parameters (with primes indicating derivatives with 
respect to the scalar field, and mp; = \/hc/G the Planck mass) as 
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which should satisfy e, |n| «& 1, the spectra can be computed using the slow-roll approximation as 

8 V 128 
A mc e A0) IA (25.5) 

3mp, € |l, jg 3mpi | pea 


In each case, the expressions on the right-hand side are to be evaluated when the scale k is equal to 
the Hubble radius during inflation. The symbol ‘~’ here indicates use of the slow-roll approximation, 
which is expected to be accurate to a few percent or better. 

From these expressions, we can compute the spectral indices [8]: 


ngc1—6e--29 ; mr —2e. (25.6) 
Another useful quantity is the ratio of the two spectra, defined by 
At) 
r= : (25.7) 
AR (Ks) 
We have 
r œ 16e c —8n, , (25.8) 


which is known as the consistency equation. 

One could consider corrections to the power-law approximation, which we discuss later. How- 
ever, for now we make the working assumption that the spectra can be approximated by such power 
laws. The consistency equation shows that r and nt are not independent parameters, and so the 
simplest inflation models give initial conditions described by three parameters, usually taken as 
AS. ns, and r, all to be evaluated at some scale k,, usually the ‘statistical center’ of the range 
explored by the data. Alternatively, one could use the parametrization V, e, and 7, all evaluated 
at a point on the putative inflationary potential. 

After the perturbations are created in the early Universe, they undergo a complex evolution 
up until the time they are observed in the present Universe. When the perturbations are small, 
this can be accurately followed using a linear theory numerical code such as CAMB or CLASS [9]. 
This works right up to the present for the CMB, but for density perturbations on small scales non- 
linear evolution is important and can be addressed by a variety of semi-analytical and numerical 
techniques. However the analysis is made, the outcome of the evolution is in principle determined 
by the cosmological model and by the parameters describing the initial perturbations, and hence 
can be used to determine them. 

Of particular interest are CMB anisotropies. Both the total intensity and two independent 
polarization modes are predicted to have anisotropies. These can be described by the radiation 
angular power spectra C; as defined in the CMB article in this volume, and again provide a complete 
description if the density perturbations are Gaussian. 
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25.1.3 The standard cosmological model 

We now have most of the ingredients in place to describe the cosmological model. Beyond those 
of the previous subsections, we need a measure of the ionization state of the Universe. The Universe 
is known to be highly ionized at redshifts below 5 or so (otherwise radiation from distant quasars 
would be heavily absorbed in the ultra-violet), and the ionized electrons can scatter microwave 
photons, altering the pattern of observed anisotropies. The most convenient parameter to describe 
this is the optical depth to scattering 7 (i.e., the probability that a given photon scatters once); in 
the approximation of instantaneous and complete reionization, this could equivalently be described 
by the redshift of reionization zi. 

As described in Sec. 25.4, models based on these parameters are able to give a good fit to 
the complete set of high-quality data available at present, and indeed some simplification is pos- 
sible. Observations are consistent with spatial flatness, and the inflation models so far described 
automatically generate negligible spatial curvature, so we can set k — 0; the density parameters 
then must sum to unity, and so one of them can be eliminated. The neutrino energy density is 
often not taken as an independent parameter; provided that the neutrino sector has the standard 
interactions, the neutrino energy density, while relativistic, can be related to the photon density 
using thermal physics arguments, and a minimal assumption takes the neutrino mass sum to be 
that of the lowest mass solution to the neutrino oscillation constraints, namely 0.06 eV. In addition, 
there is no observational evidence for the existence of tensor perturbations (with the upper limits 
now starting to become constraining on models), and so r could be set to zero. This leaves seven 
parameters, which is the smallest set that can usefully be compared to the present cosmological 
data. This model is referred to by various names, including ACDM, the concordance cosmology, 
and the standard cosmological model. 

Of these parameters, only Q, is accurately measured directly. The radiation density is dom- 
inated by the energy in the CMB, and the COBE satellite FIRAS experiment determined its 
temperature to be T = 2.7255 + 0.0006 K [10],! corresponding to Q, = 2.47 x 10 ?A ?. It typi- 
cally can be taken as fixed when fitting other data. Hence the minimum number of cosmological 
parameters varied in fits to data is six, though as described below there may additionally be many 
‘nuisance’ parameters necessary to describe astrophysical processes influencing the data. 

In addition to this minimal set, there is a range of other parameters that might prove important 
in future as the data-sets further improve, but for which there is so far no direct evidence, allowing 
them to be set to specific values for now. We discuss various speculative options in the next 
section. For completeness at this point, we mention one other interesting quantity, the helium 
fraction, which is a non-zero parameter that can affect the CMB anisotropies at a subtle level. It is 
usually fixed in microwave anisotropy studies, but the data are approaching a level where allowing 
its variation may become mandatory. 

In conventional parameter estimation, a set of parameters is chosen by hand and the aim is 
to constrain their values. The higher-level inference problem of model selection instead compares 
different choices of parameter sets, as is necessary to assess whether observations are pointing 
towards inclusion of new physical effects. Bayesian inference offers an attractive framework for 
cosmological model selection, setting a tension between model predictiveness and ability to fit the 
data [11], and its use is becoming widespread. 


25.1.4 Derived parameters 
'The parameter list of the previous subsection is sufficient to give a complete description of cos- 
mological models that agree with observational data. However, it is not a unique parameterization, 
‘All quoted uncertainties in this article are 10/6896 confidence and all upper limits are 95% confidence. Cosmolog- 


ical parameters sometimes have significantly non-Gaussian uncertainties. Throughout we have rounded central values, 
and especially uncertainties, from original sources, in cases where they appear to be given to excessive precision. 
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and one could instead use parameters derived from that basic set. Parameters that can be obtained 
from the set given above include the age of the Universe, the present horizon distance, the present 
neutrino background temperature, the epoch of matter-radiation equality, the epochs of recombi- 
nation and decoupling, the epoch of transition to an accelerating Universe, the baryon-to-photon 
ratio, and the baryon-to-dark-matter density ratio. In addition, the physical densities of the matter 
components, Q;h?, are often more useful than the density parameters. The density perturbation 
amplitude can be specified in many different ways other than the large-scale primordial amplitude, 
for instance, in terms of its effect on the CMB, or by specifying a short-scale quantity, a common 
choice being the present linear-theory mass dispersion on a radius of 8 h-! Mpc, known as og. 

Different types of observation are sensitive to different subsets of the full cosmological parameter 
set, and some are more naturally interpreted in terms of some of the derived parameters of this 
subsection than on the original base parameter set. In particular, most types of observation feature 
degeneracies whereby they are unable to separate the effects of simultaneously varying specific 
combinations of several of the base parameters. 


25.2 Extensions to the standard model 

At present, there is no positive evidence in favor of extensions of the standard model. These 
are becoming increasingly constrained by the data, though there always remains the possibility of 
trace effects at a level below present observational capability. 


25.2.1 More general perturbations 

The standard cosmology assumes adiabatic, Gaussian perturbations. Adiabaticity means that 
all types of material in the Universe share a common perturbation, so that if the space-time is 
foliated by constant-density hypersurfaces, then all fluids and fields are homogeneous on those slices, 
with the perturbations completely described by the variation of the spatial curvature of the slices. 
Gaussianity means that the initial perturbations obey Gaussian statistics, with the amplitudes 
of waves of different wavenumbers being randomly drawn from a Gaussian distribution of width 
given by the power spectrum. Note that gravitational instability generates non-Gaussianity; in this 
context, Gaussianity refers to a property of the initial perturbations, before they evolve. 

The simplest inflation models, based on one dynamical field, predict adiabatic perturbations 
and a level of non-Gaussianity that is too small to be detected by any experiment so far conceived. 
For present data, the primordial spectra are usually assumed to be power laws. 


25.2.1.1 Non-power-law spectra 

For typical inflation models, it is an approximation to take the spectra as power laws, albeit 
usually a good one. As data quality improves, one might expect this approximation to come under 
pressure, requiring a more accurate description of the initial spectra, particularly for the density 
perturbations. In general, one can expand In AL as 


kl dne| Lok 
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In AZ (k) = In A2 (kx) + (ns. — 1) In Te, (25.9) 
where the coefficients are all evaluated at some scale ką. The term dn;/dlnk|, is often called 
the running of the spectral index [12]. Once non-power-law spectra are allowed, it is necessary to 
specify the scale k, at which the spectral index is defined. 


25.2.1.2 Isocurvature perturbations 

An isocurvature perturbation is one that leaves the total density unperturbed, while perturbing 
the relative amounts of different materials. If the Universe contains N fluids, there is one growing 
adiabatic mode and N — 1 growing isocurvature modes (for reviews see Ref. [7] and Ref. [13]). 
'These can be excited, for example, in inflationary models where there are two or more fields that 
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acquire dynamically-important perturbations. If one field decays to form normal matter, while 
the second survives to become the dark matter, this will generate a cold dark matter isocurvature 
perturbation. 

In general, there are also correlations between the different modes, and so the full set of per- 
turbations is described by a matrix giving the spectra and their correlations. Constraining such a 
general construct is challenging, though constraints on individual modes are beginning to become 
meaningful, with no evidence that any other than the adiabatic mode must be non-zero. 


25.2.1.3 Seeded perturbations 

An alternative to laying down perturbations at very early epochs is that they are seeded through- 
out cosmic history, for instance by topological defects such as cosmic strings. It has long been 
excluded that these are the sole original of structure, but they could contribute part of the pertur- 
bation signal, current limits being just a few percent [14]. In particular, cosmic defects formed in 
a phase transition ending inflation is a plausible scenario for such a contribution. 


25.2.1.4 Non- Gaussianity 

Multi-field inflation models can also generate primordial non-Gaussianity (reviewed, e.g., in 
Ref. |7]). The extra fields can either be in the same sector of the underlying theory as the inflaton, 
or completely separate, an interesting example of the latter being the curvaton model [15]. Current 
upper limits on non-Gaussianity are becoming stringent, but there remains strong motivation to 
push down those limits and perhaps reveal trace non-Gaussianity in the data. If non-Gaussianity is 
observed, its nature may favor an inflationary origin, or a different one such as topological defects. 


25.2.2 Dark matter properties 

Dark matter properties are discussed in the Dark Matter chapter in this volume. The simplest 
assumption concerning the dark matter is that it has no significant interactions with other matter, 
and that its particles have a negligible velocity as far as structure formation is concerned. Such 
dark matter is described as ‘cold,’ and candidates include the lightest supersymmetric particle, the 
axion, and primordial black holes. As far as astrophysicists are concerned, a complete specification 
of the relevant cold dark matter properties is given by the density parameter Qe, though those 
seeking to detect it directly need also to know its interaction properties. 

Cold dark matter is the standard assumption and gives an excellent fit to observations, except 
possibly on the shortest scales where there remains some controversy concerning the structure of 
dwarf galaxies and possible substructure in galaxy halos. It has long been excluded for all the dark 
matter to have a large velocity dispersion, so-called ‘hot’ dark matter, as it does not permit galaxies 
to form; for thermal relics the mass must be above about 1 keV to satisfy this constraint, though 
relics produced non-thermally, such as the axion, need not obey this limit. However, in future 
further parameters might need to be introduced to describe dark matter properties relevant to 
astrophysical observations. Suggestions that have been made include a modest velocity dispersion 
(warm dark matter) and dark matter self-interactions. There remains the possibility that the dark 
matter is comprized of two separate components, e.g., a cold one and a hot one, an example being 
if massive neutrinos have a non-negligible effect. 


25.2.3 Relativistic species 

The number of relativistic species in the young Universe (omitting photons) is denoted Neg. In 
the standard cosmological model only the three neutrino species contribute, and its baseline value 
is assumed fixed at 3.044 (the small shift from 3 is because of a slight predicted deviation from a 
thermal distribution [16]). However other species could contribute, for example an extra neutrino, 
possibly of sterile type, or massless Goldstone bosons or other scalars. It is hence interesting to 
study the effect of allowing this parameter to vary, and indeed although 3.044 is consistent with 
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the data, most analyses currently suggest a somewhat higher value (e.g., Ref. [17]). 


25.2.4 Dark energy and modified gravity 

While the standard cosmological model given above features a cosmological constant, in order to 
explain observations indicating that the Universe is presently accelerating, further possibilities exist 
under the general headings of ‘dark energy’ and ‘modified gravity. These topics are described in 
detail in the Dark Energy chapter in this volume. This article focuses on the case of the cosmological 
constant, since this simple model is a good match to existing data. We note that more general 
treatments of dark energy/modified gravity will lead to weaker constraints on other parameters. 


25.2.5 Complex ionization history 

The full ionization history of the Universe is given by the ionization fraction as a function 
of redshift z. The simplest scenario takes the ionization to have the small residual value left 
after recombination up to some redshift z;, at which point the Universe instantaneously reionizes 
completely. Then there is a one-to-one correspondence between 7 and z (that relation, however, 
also depending on other cosmological parameters). An accurate treatment of this process will track 
separate histories for hydrogen and helium. While currently rapid ionization appears to be a good 
approximation, as data improve a more complex ionization history may need to be considered. 


25.2.6 Varying ‘constants’ 

Variation of the fundamental constants of Nature over cosmological times is another possible 
enhancement of the standard cosmology. There is a long history of study of variation of the 
gravitational constant Gy, and more recently attention has been drawn to the possibility of small 
fractional variations in the fine-structure constant. There is presently no observational evidence for 
the former, which is tightly constrained by a variety of measurements. Evidence for the latter has 
been claimed from studies of spectral line shifts in quasar spectra at redshift z ~ 2 [18], but this is 
presently controversial and in need of further observational study. 


25.2.7 Cosmic topology 

The usual hypothesis is that the Universe has the simplest topology consistent with its geometry, 
for example that a flat Universe extends forever. Observations cannot tell us whether that is true, 
but they can test the possibility of a non-trivial topology on scales up to roughly the present Hubble 
scale. Extra parameters would be needed to specify both the type and scale of the topology; for 
example, a cuboidal topology would need specification of the three principal axis lengths and 
orientation. At present, there is no evidence for non-trivial cosmic topology [19]. 


25.3 Cosmological Probes 

The goal of the observational cosmologist is to utilize astronomical information to derive cos- 
mological parameters. The transformation from the observables to the parameters usually involves 
many assumptions about the nature of the data, as well as of the dark sector. Below we outline the 
physical processes involved in each of the major probes, and the main recent results. The first two 
subsections concern probes of the homogeneous Universe, while the remainder consider constraints 
from perturbations. 

In addition to statistical uncertainties we note three sources of systematic uncertainties that will 
apply to the cosmological parameters of interest: (i) due to the assumptions on the cosmological 
model and its priors (¢.e., the number of assumed cosmological parameters and their allowed range); 
(ii) due to the uncertainty in the astrophysics of the objects (e.g., light-curve fitting for supernovae 
or the mass-temperature relation of galaxy clusters); and (iii) due to instrumental and observational 
limitations (e.g., the effect of ‘seeing’ on weak gravitational lensing measurements, or beam shape 
on CMB anisotropy measurements). 

These systematics, the last two of which appear as ‘nuisance parameters’, pose a challenging 
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problem to the statistical analysis. We attempt a statistical fit to the whole Universe with 6 
to 12 parameters, but we might need to include hundreds of nuisance parameters, some of them 
highly correlated with the cosmological parameters of interest (for example time-dependent galaxy 
biasing could mimic the growth of mass fluctuations). Fortunately, there is some astrophysical 
prior knowledge on these effects, and a small number of physically-motivated free parameters would 
ideally be preferred in the cosmological parameter analysis. 


25.3.1 Measures of the Hubble constant 

In 1929, Edwin Hubble discovered the law of expansion of the Universe by measuring dis- 
tances to nearby galaxies. The slope of the relation between the distance and recession veloc- 
ity is defined to be the present-epoch Hubble constant, Ho. Astronomers argued for decades 
about the systematic uncertainties in various methods and derived values over the wide range 
40kms-! Mpc^! $ Ho <£ 100kms~! Mpc. 

One of the most reliable results on the Hubble constant came from the Hubble Space Telescope 
(HST) Key Project [20]. This study used the empirical period-luminosity relation for Cepheid 
variable stars, and calibrated a number of secondary distance indicators: Type Ia Supernovae 
(SNe Ia), the Tully—Fisher relation, surface-brightness fluctuations, and Type II Supernovae. This 
approach has been further extended. Based on HST photometry of 75 Milky Way Cepheids and 
Gaia EDR3 parallaxes, the SHOES team derived Hy = 73.2 + 1.3 kms~! Mpc™? [21]. 

Three other methods have been used recently. One is a calibration of the tip of the red-giant 
branch applied to Type Ia supernovae, the Carnegie-Chicago Hubble Programme (CCHP) finding 
Ho = 69.8 + 0.6 (stat.) + 1.6 (sys. kms^!Mpc ! [22]. The second uses the method of time 
delay in gravitationally-lensed quasars; Birrer et al. [23] find Hp = 67.4133 kms-! Mpc^! from a 
sample of 40 lenses (which includes as a subset the six systems previously used by the HOLICOW 
collaboration [24] who obtained a higher value). A third method that came to fruition recently 
is based on gravitational waves; the ‘bright standard siren’ method applied to the binary neutron 
star GW170817 yields Ho = "ys kms-!Mpc ! [25]. Adding two ‘dark standard siren’ systems 
shifts this to Hy = a km s~! Mpc^! [26], still dominated by the single bright siren system. 
When many more gravitational-wave events have been acquired, the future uncertainties on Ho 
from standard sirens will get smaller. 


The determination of Hp by the Planck Collaboration [2] gives a lower value than most of the 
above methods, Hp = 67.4 + 0.5 kms^! Mpc^!. As they discuss, there is strong degeneracy of Ho 
with other parameters, particularly Qm and the neutrino mass. It is worth noting that using the 
‘inverse distance ladder’ method gives a result Hy = 67.8 + 1.3 km s^! Mpc ! [27], close to the 
Planck result. The inverse distance ladder relies on absolute-distance measurements from baryon 
acoustic oscillations (BAOs) to calibrate the intrinsic magnitude of the SNe Ia (rather than nearby 
Cepheids and parallax). This measurement was derived from 207 spectroscopically-confirmed Type 
Ia supernovae from the Dark Energy Survey (DES), an additional 122 low-redshift SNe Ia, and 
measurements of BAOs. A combination of DES Year 3 (Y3) clustering and weak lensing with BAO 
and BBN (assuming ACDM) gives Ho = 67.6 + 0.9 kms~! Mpc! [28]. The completed Extended 
Baryon Oscillation Spectroscopic Survey (eBOSS) [29] inverse distance ladder result, within an 
assumed extended cosmological model, is Hp = 68.2 + 0.8 km s^! Mpc ^ !, also close to the Planck 
value. 


The tension between the Hp values from Planck and the traditional cosmic distance ladder 
methods, with the SHOES result deviating from Planck by 4.20, is under intense investigation 
for potential systematic effects. There is possibly a trend for higher Hop derived from the nearby 
Universe and a lower Ho from the early Universe, which has led some researchers to propose a time- 
variation of the dark energy component or other exotic scenarios. Ongoing studies are addressing 
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the question of whether the Hubble tension is due to systematics in at least one of the probes, or 
a signature of new physics. 

Figure 25.1 shows a selection of recent Ho values, summarizing the current status of the Hub- 
ble constant tension. We note that while the tension remains highly significant, its severity has 
somewhat lessened compared to our assessment in the previous edition of this article two years ago. 
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Figure 25.1: A selection of recent Hp measurements from the various projects as described in 
the text. The standard-siren determinations are omitted as they are too wide for the plot. Figure 
courtesy of Pablo Lemos. 


25.3.2 Supernovae as cosmological probes 

Empirically, the peak luminosity of SNe Ia can be used as an efficient distance indicator 
(e.g., Ref. [30]), thus allowing cosmology to be constrained via the distance-redshift relation. The 
favorite theoretical explanation for SNe Ia is the thermonuclear disruption of carbon-oxygen white 
dwarfs. Although not perfect ‘standard candles’, it has been demonstrated that by correcting for 
a relation between the light-curve shape, color, and luminosity at maximum brightness, the dis- 
persion of the measured luminosities can be greatly reduced. T'here are several possible systematic 
effects that may affect the accuracy of the use of SNe Ia as distance indicators, e.g., evolution with 
redshift and interstellar extinction in the host galaxy and in the Milky Way. 

'Two major studies, the Supernova Cosmology Project and the High-z Supernova Search Team, 
found evidence for an accelerating Universe [31], interpreted as due to a cosmological constant or 
a dark energy component. When combined with the CMB data (which indicate near flatness, i.e., 
Om + O4 ~ 1), the best-fit values were Om œ 0.3 and Q4 œ 0.7. Most results in the literature 
are consistent with the w = —1 cosmological constant case. The most-used sample currently 
is the Pantheon compilation [32]. This sample of 1048 spectroscopically-confirmed SNe Ia gives 
Om = 0.298 + 0.022 (stat+sym) for an assumed flat ACDM model. In combination with the 
CMB, for a flat wCDM model these data give w = —1.03 + 0.04 and Êm = 0.307 + 0.012. For 
comparison, an analysis of a sample of 207 spectroscopically-confirmed DES SNe Ia combined with 
122 low-redshift SNe [33] yielded Qm = 0.331 + 0.038 for an assumed flat ACDM model. Future 
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experiments will refine constraints on the cosmic equation of state w(z). 


25.3.3 Cosmic microwave background 

The physics of the CMB is described in detail in the CMB chapter in this volume. Before 
recombination, the baryons and photons are tightly coupled, and the perturbations oscillate in 
the potential wells generated primarily by the dark matter perturbations. After decoupling, the 
baryons are free to collapse into those potential wells. The CMB carries a record of conditions at 
the time of last scattering, often called primary anisotropies. In addition, it is affected by various 
processes as it propagates towards us, including the effect of a time-varying gravitational potential 
(the integrated Sachs-Wolfe effect), gravitational lensing, and scattering from ionized gas at low 
redshift. 

The primary anisotropies, the integrated Sachs-Wolfe effect, and the scattering from a homoge- 
neous distribution of ionized gas, can all be calculated using linear perturbation theory. Available 
codes include CAMB and CLASS [9], the former widely used embedded within the analysis package 
CosmoMC [34] and in higher-level analysis packages such as CosmoSIS [35] and CosmoLike [36]. Grav- 
itational lensing is also calculated in these codes. Secondary effects, such as inhomogeneities in the 
reionization process, and scattering from gravitationally-collapsed gas (the Sunyaev-Zeldovich or 
SZ effect), require more complicated, and more uncertain, calculations. 

The upshot is that the detailed pattern of anisotropies depends on all of the cosmological 
parameters. In a typical cosmology, the anisotropy power spectrum [usually plotted as (£ + 1)C;] 
features a flat plateau at large angular scales (small £), followed by a series of oscillatory features 
at higher angular scales, the first and most prominent being at around one degree (/ ~ 200). These 
features, known as acoustic peaks, represent the oscillations of the photon-baryon fluid around 
the time of decoupling. Some features can be closely related to specific parameters—for instance, 
the location in multipole space of the set of peaks probes the spatial geometry, while the relative 
heights of the peaks probe the baryon density—but many other parameters combine to determine 
the overall shape. 

The 2018 data release from the Planck satellite [1] gives the most powerful results to date on 
the spectrum of CMB temperature anisotropies, with a precision determination of the temperature 
power spectrum to beyond £ = 2000. The Atacama Cosmology Telescope (ACT) and South Pole 
Telescope (SPT) experiments extend these results to higher angular resolution, though without 
full-sky coverage. Planck and the polarization-sensitive versions of ACT and SPT give the state 
of the art in measuring the spectrum of E-polarization anisotropies and the correlation spectrum 
between temperature and polarization. These are consistent with models based on the parameters 
we have described, and provide accurate determinations of many of those parameters [2]. Primordial 
B-mode polarization has not been detected (although the gravitational lensing effect on B modes 
has been measured). 

The data provide an exquisite measurement of the location of the set of acoustic peaks, de- 
termining the angular-diameter distance of the last-scattering surface. In combination with other 
data this strongly constrains the spatial geometry, in a manner consistent with spatial flatness and 
excluding significantly-curved Universes. CMB data give a precision measurement of the age of the 
Universe. The CMB also gives a baryon density consistent with, and at higher precision than, that 
coming from BBN. It affirms the need for both dark matter and dark energy. It shows no evidence 
for dynamics of the dark energy, being consistent with a pure cosmological constant (w = —1). 
'The density perturbations are consistent with a power-law primordial spectrum, and there is no 
indication yet of tensor perturbations. The current best-fit for the reionization optical depth from 
CMB data, 7 = 0.054, is in line with models of how early structure formation induces reionization. 

Planck has also made the first all-sky map of the CMB lensing field, which probes the entire 
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matter distribution in the Universe and adds some additional constraining power to the CMB-only 
data-sets. These measurements are compatible with the expected effect in the standard cosmology. 


25.3.4 Galaxy clustering 

The power spectrum of density perturbations is affected by the nature of the dark matter. 
Within the ACDM model, the power spectrum shape depends primarily on the primordial power 
spectrum and on the combination mAh, which determines the horizon scale at matter—radiation 
equality, with a subdominant dependence on the baryon density. The matter distribution is most 
easily probed by observing the galaxy distribution, but this must be done with care since the galaxies 
do not perfectly trace the dark matter distribution. Rather, they are a ‘biased’ tracer of the dark 
matter [37]. The need to allow for such bias is emphasized by the observation that different types of 
galaxies show bias with respect to each other. In particular, scale-dependent and stochastic biasing 
may introduce a systematic effect on the determination of cosmological parameters from redshift 
surveys [38]. Prior knowledge from simulations of galaxy formation or from gravitational lensing 
data could help to quantify biasing. Furthermore, the observed 3D galaxy distribution is in redshift 
space, i.e., the observed redshift is the sum of the Hubble expansion and the line-of-sight peculiar 
velocity, leading to linear and non-linear dynamical effects that also depend on the cosmological 
parameters. On the largest length scales, the galaxies are expected to trace the location of the 
dark matter, except for a constant multiplier b to the power spectrum, known as the linear bias 
parameter. On scales smaller than 20 Mpc or so, the clustering pattern is ‘squashed’ in the radial 
direction due to coherent infall, which depends approximately on the parameter 8 = 0°-°/b (on 
these shorter scales, more complicated forms of biasing are not excluded by the data). On scales of 
a few Mpc, there is an effect of elongation along the line of sight (colloquially known as the 'finger 
of God' effect) that depends on the galaxy velocity dispersion. 


25.3.4.1 Baryon acoustic oscillations 

The power spectra of the 2-degree Field (2dF) Galaxy Redshift Survey and the Sloan Digital 
Sky Survey (SDSS) are well fit by a ACDM model and both surveys showed first evidence for baryon 
acoustic oscillations (BAOs) [39,40]. When eBOSS is combined with Planck, Pantheon Type Ia 
Supernovae and other probes the result is wp = —1.020 + 0.032 at the pivot redshift zp = 0.29 [29]. 
Similar results for w were obtained e.g. by the WiggleZ survey [41]. 


25.3.4.2 Redshift distortion 

There is continuing interest in the ‘redshift distortion’ effect. This distortion depends on cos- 
mological parameters [42] via the perturbation growth rate in linear theory f(z) = dlInó/dlna ~ 
QX (z), where y ~ 0.55 for the ACDM model and may be different for modified gravity models. By 
measuring f(z) it is feasible to constrain y and rule out certain modified gravity models [43,44]. We 
note the degeneracy of the redshift-distortion pattern and the geometric distortion (the so-called 
Alcock-Paczynski effect [45]), e.g., as illustrated by the WiggleZ survey [46] and eBOSS [29]. 


25.3.4.3 Limits on neutrino mass from galaxy surveys and other probes 

Large-scale structure data place constraints on Q, due to the neutrino free-streaming effect [47]. 
Presently there is no clear detection, and upper limits on neutrino mass are commonly estimated 
by comparing the observed galaxy power spectrum with a four-component model of baryons, cold 
dark matter, a cosmological constant, and massive neutrinos. Such analyses also assume that the 
primordial power spectrum is adiabatic, scale-invariant, and Gaussian. Potential systematic effects 
include biasing of the galaxy distribution and non-linearities of the power spectrum. An upper limit 
can also be derived from CMB anisotropies alone, while combination with additional cosmological 
data-sets can improve the results. 

'The most recent results on neutrino mass upper limits and other neutrino properties are summa- 
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rized in the Neutrinos in Cosmology chapter in this volume. The latest cosmological data constrain 
the sum of neutrino masses to be below 0.2 eV or even 0.1 eV depending on the assumed priors 
and systematics (e.g. Refs. [2,29].) Since the lower limit on this sum from oscillation experiments 
(assuming normal hierarchy) is 0.06 eV it is expected that future cosmological surveys will soon 
detect effects from the neutrino mass. Also, current cosmological datasets are in good agreement 
with the standard value for the effective number of neutrino species Neg = 3.044. 


25.3.5 Clustering in the inter-galactic medium 


It is commonly assumed, based on hydrodynamic simulations, that the neutral hydrogen in the 
inter-galactic medium (IGM) can be related to the underlying mass distribution. It is then possible 
to estimate the matter power spectrum from the absorption observed in quasar spectra, the so-called 
Lyman-o forest. The usual procedure is to measure the power spectrum of the transmitted flux, 
and then to infer the mass power spectrum. Photo-ionization heating by the ultraviolet background 
radiation and adiabatic cooling by the expansion of the Universe combine to give a simple power- 
law relation between the gas temperature and the baryon density. It also follows that there is a 
power-law relation between the optical depth 7 and pp. Therefore, the observed flux F = exp(—7) 
is strongly correlated with pp, which itself traces the mass density. The matter and flux power 
spectra can be related by a biasing function that is calibrated from simulations. There are two 
variants of Lyman-alpha analyses: 1-dimensional power spectra from individual lines-of-sight that 
probe small (~Mpc) scales, and 3-dimensional Lyman-alpha BAO analyses that measure large-scale 
correlations (over ~100 Mpc scales) using neighbouring quasar lines-of-sight. 


The latest BOSS and eBOSS datasets have provided measurements of the BAO scale both in 
the Lyman-q@ absorption and in its cross correlation with quasars at an effective redshift z = 2.3 
[29,48]. The results agree within 1.5-o with the Planck CMB flat-A estimation for BAO scale. 
Such measurements will improve with the Dark Energy Spectroscopic Instrument (DESI) and other 
future spectroscopic surveys. The Lyman-o flux power spectrum has also been used to constrain 
the nature of dark matter, for example limiting the amount of warm dark matter [49]. 


25.3.6 Weak gravitational lensing 


Images of background galaxies are distorted by the gravitational effect of mass variations along 
the line of sight. Deep gravitational potential wells, such as galaxy clusters, generate ‘strong lensing’ 
leading to arcs, arclets, and multiple images, while more moderate perturbations give rise to ^weak 
lensing’, Weak lensing is now widely used to measure the mass power spectrum in selected regions 
of the sky (see Ref. [50] for reviews). Since the signal is weak, the image of deformed galaxy shapes 
(the ‘shear map’) must be analyzed statistically to measure the power spectrum, higher moments, 
and cosmological parameters. There are various systematic effects in the interpretation of weak 
lensing, e.g., due to atmospheric distortions during observations, the redshift distribution of the 
background galaxies (usually depending on the accuracy of photometric redshifts), the intrinsic 
correlation of galaxy shapes, and non-linear modeling uncertainties. 


Weak-lensing measurements from the Kilo-Degree Survey (KiDS) [51], the Subaru Hyper- 
Suprime-Cam (HSC) [52], and DES Y3 [53] have constrained the clumpiness parameter Sg = 
o8(Qm/0.3)°°. Each of these surveys has yielded Sg lower by about 2c than Sg derived from 
Planck. This tension is not yet resolved. Figure 25.2 shows the €4,—5& constraints derived from 
of KiDS (weak lensing) and DES Y3 (weak lensing with and without galaxy clustering) versus the 
CMB constraint from Planck. Results from weak lensing from DES, combined with other probes, 
are described in the next section. 


24th January, 2022 


13 25. Cosmological Parameters 


0.900 


DES 3x2pt 
DES + 

KiDS COSEBIs 
Planck CMB 


0.875 
0.850 
0.825 


0.800 


S 


0.775 


0.750 


0.725 


0.700 


ders 0.1 0.2 0.3 0.4 0.5 


Figure 25.2: Marginalized posterior contours (inner 68% confidence level, outer 95% confidence 
level) in the Q,,—$s plane. The plot shows a subset of the cases shown in Figure 19 of Ref. [28]: 
KiDS weak lensing alone, DES Y3 weak lensing alone, DES Y3 combined weak lensing and galaxy 
clustering (3x2pt), and Planck CMB. Figure courtesy of Michael Troxel and the DES collaboration. 


25.3.7 Other probes 

Other probes that have been used to constrain cosmological parameters, but that are not 
presently competitive in terms of accuracy, are the integrated Sachs-Wolfe effect [54,55], the num- 
ber density or composition of galaxy clusters [56], and galaxy peculiar velocities which probe the 
mass fluctuations in the local Universe [57]. 


25.4 Bringing probes together 

Although it contains two ingredients—dark matter and dark energy—which have not yet been 
verified by laboratory experiments, the ACDM model is almost universally accepted by cosmologists 
as the best description of the present data. The approximate values of some of the key parameters 
are Q œ 0.05, Q.e œ 0.25, Qa œ 0.70, and a Hubble constant h œ 0.70. The spatial geometry is 
very close to flat (and usually assumed to be precisely flat), and the initial perturbations Gaussian, 
adiabatic, and nearly scale-invariant. 

'The most powerful data source is the CMB, which on its own supports all these main tenets. 
Values for some parameters, as given in Ref. [2], are reproduced in Table 25.1. These particular 
results presume a flat Universe. The constraints are somewhat strengthened by adding additional 
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data-sets, BAO being shown in the Table as an example, though most of the constraining power 
resides in the CMB data. Similar constraints at lower precision were previously obtained by the 
WMAP collaboration. 


Table 25.1: Parameter constraints reproduced from Ref. [2] (Table 2, 
column 5), with some additional rounding. Both columns assume the 
ACDM cosmology with a power-law initial spectrum, no tensors, spatial 
flatness, a cosmological constant as dark energy, and the sum of neutrino 
masses fixed to 0.06 eV. Above the line are the six parameter combinations 
actually fit to the data (Ac is a measure of the sound horizon at last 
scattering); those below the line are derived from these. The first column 
uses Planck primary CMB data plus the Planck measurement of CMB 
lensing. This column gives our present recommended values. The second 
column adds in data from a compilation of BAO measurements described 
in Ref. [2]. The perturbation amplitude A2, (denoted A, in the original 
paper) is specified at the scale 0.05 Mpc !. Uncertainties are shown at 
6896 confidence. 


Planck TT,TE,EE+lowE+lensing +BAO 


Qy? 0.02237 + 0.00015 0.02242 + 0.00014 
Qh? 0.1200 + 0.0012 0.1193 + 0.0009 
100 Qc 1.0409 + 0.0003 1.0410 + 0.0003 
Ns 0.965 + 0.004 0.966 + 0.004 
T 0.054 + 0.007 0.056 + 0.007 
In(1019 A2) 3.044 + 0.014 3.047 + 0.014 
h 0.674 + 0.005 0.677 + 0.004 
0g 0.811 + 0.006 0.810 + 0.006 
Ou 0.315 + 0.007 0.311 + 0.006 
QA 0.685 + 0.007 0.689 + 0.006 


If the assumption of spatial flatness is lifted, it turns out that the primary CMB on its own con- 
strains the spatial curvature fairly weakly, due to a parameter degeneracy in the angular-diameter 
distance. However, inclusion of other data readily removes this degeneracy. Simply adding the 
Planck lensing measurement, and with the assumption that the dark energy is a cosmological 
constant, yields a 68% confidence constraint on Oot = X Q; + Qa = 1.011 x 0.006 and further 
adding BAO makes it 0.9993 + 0.0019 [2]. Results of this type are normally taken as justifying the 
restriction to flat cosmologies. 


One derived parameter that is very robust is the age of the Universe, since there is a useful 
coincidence that for a flat Universe the position of the first peak is strongly correlated with the 
age. The CMB data give 13.797 + 0.023 Gyr (assuming flatness). This is in good agreement with 
the ages of the oldest globular clusters and with radioactive dating. 

The baryon density Qp is now measured with high accuracy from CMB data alone, and is 
consistent with and more precise than the determination from BBN. The value quoted in the 
Big-Bang Nucleosynthesis chapter in this volume is QA? = 0.0224 + 0.0007. 

While Q, is measured to be non-zero with very high confidence, there is no evidence of evolution 
of the dark energy density. As described in the Dark Energy chapter in this volume a combination 
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of CMB, SN, and BAO measurements, assuming a flat Universe, found w = —1.03 + 0.03 [29], 
consistent with the cosmological constant case w = —1. Allowing more complicated forms of dark 
energy weakens the limits. 


The data provide strong support for the main predictions of the simplest inflation models: 
spatial flatness and adiabatic, Gaussian, nearly scale-invariant density perturbations. But it is 
disappointing that there is no sign of primordial gravitational waves; combining Planck and WMAP 
with BICEP2/Keck Array BK18 data (plus BAO data to help constrain ns) gives a 95% confidence 
upper limit of r < 0.036 at the scale 0.05 Mpc ^! [58]. The density perturbation spectral index is 
clearly required to be less than one by current data, though the strength of that conclusion can 
weaken if additional parameters are included in the model fits. 

Tests have been made for various types of non-Gaussianity, a particular example being a param- 
eter fwr, that measures a quadratic contribution to the perturbations. Various non-Gaussian shapes 
are possible (see Ref. [59] for details), and current constraints on the popular ‘local’, ‘equilateral’, 


and ‘orthogonal’ types (combining CMB temperature and polarization data) are local = —-1+5, 
D = —26 + 47, and fa = —38 + 24 respectively (these look weak, but prominent non- 


Gaussianity requires the product fNLAR to be large, and Ap is of order 1075). Clearly none of 
these give any indication of primordial non-Gaussianity. 

While the above results come from the CMB alone, other probes are becoming competitive 
(especially when considering more complex cosmological models), and so combination of data from 
different sources is of growing importance. We note that it has become fashionable to combine 
probes at the level of power-spectrum data vectors, taking into account nuisance parameters in 
each type of measurement. Discussions on ‘tension’ in resulting cosmological parameters depend 
on the statistical approaches used. Commonly the cosmology community works within the Bayesian 
framework, and assesses agreement amongst data sets with respect to a model via Bayesian Evi- 
dence, essentially the denominator in Bayes's theorem. 

As an example of results, combining DES Y3 (position-position clustering, galaxy-galaxy lens- 
ing, and weak lensing shear) with Planck, BAO and RSD measurements from eBOSS, and type Ia su- 
pernovae from the Pantheon dataset compilation has shown the datasets to be mutually compatible 
and yields very tight constraints on cosmological parameters: Sg = og(Qm/0.3)°> = 0.812 + 0.008, 
and Qn = 0.306790}; in ACDM, and w = —1.03 + 0.03 in wCDM [28] matching the con- 
straint in Ref. [29]. The combined measurement of the Hubble constant within ACDM gives 
Ho = 65.0275 kms! Mpc 1, still leaving substantial tension with the SHOES measurement de- 
scribed earlier. Future analyses and the next generation of surveys will test for deviations from 
ACDM, for example epoch-dependent w(z) and modifications to General Relativity. 


25.5 Outlook for the future 


The concordance model is now well-established, and there seems little room left for any dramatic 
revision of this paradigm. A measure of the strength of that statement is how difficult it has proven 
to formulate convincing alternatives. 


Should there indeed be no major revision of the current paradigm, we can expect future de- 
velopments to take one of two directions. Either the existing parameter set will continue to prove 
sufficient to explain the data, with the parameters subject to ever-tightening constraints, or it 
will become necessary to deploy new parameters. The latter outcome would be very much the 
more interesting, offering a route towards understanding new physical processes relevant to the 
cosmological evolution. There are many possibilities on offer for striking discoveries, for example: 


e the cosmological effects of a neutrino mass may be unambiguously detected, shedding light 
on fundamental neutrino properties; 
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e detection of primordial non-Gaussianities would indicate that non-linear processes influence 
the perturbation generation mechanism; 

e detection of variation in the dark-energy density (i.e., w Æ —1) would provide much-needed 
experimental input into its nature. 


These provide more than enough motivation for continued efforts to test the cosmological model 
and improve its accuracy. Over the coming years, there are a wide range of new observations that 
will bring further precision to cosmological studies. Indeed, there are far too many for us to be 
able to mention them all here, and so we will just highlight a few areas. 

The CMB observations will improve in several directions. A current frontier is the study of 
polarization, for which power spectrum measurements have now been made by several experiments. 
Detection of primordial B-mode anisotropies is the next major goal and a variety of projects are 
targeting this, though theory gives little guidance as to the likely signal level. Future CMB projects 
that are approved include LiteBIRD and the Simons Observatory. 

An impressive array of cosmology surveys are already operational, under construction, or pro- 
posed, including the ground-based Hyper Suprime Camera (HSC) and Rubin-LSST imaging sur- 
veys, spectroscopic surveys such as DESI and space missions Euclid and the Roman Wide-Field 
Infrared Survey (WFIRST). 

An exciting area for the future is radio surveys of the redshifted 21-cm line of hydrogen. Be- 
cause of the intrinsic narrowness of this line, by tuning the bandpass the emission from narrow 
redshift slices of the Universe will be measured to extremely high redshift, probing the details of 
the reionization process at redshifts up to perhaps 20, as well as measuring large-scale features 
such as the BAOs. LOFAR and CHIME are the first instruments able to do this and have begun 
operations. In the longer term, the Square Kilometre Array (SKA) will take these studies to a 
precision level. 

The development of the first precision cosmological model is a major achievement. However, it is 
important not to lose sight of the motivation for developing such a model, which is to understand the 
underlying physical processes at work governing the Universe's evolution. From that perspective, 
progress has been much less dramatic. For instance, there are many proposals for the nature of 
the dark matter, but no consensus as to which is correct. The nature of the dark energy remains a 
mystery. Even the baryon density, now measured to an accuracy of a percent, lacks an underlying 
theory able to predict it within orders of magnitude. Precision cosmology may have arrived, but 
at present many key questions remain to motivate and challenge the cosmology community. 
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